Multiplicity of X-Ray Selected T Tauri Stars in Chamaeleon 
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ABSTRACT 

We report on a multiplicity survey of a sample of X-ray selected young stars in the Cha- 
maeleon association. We used speckle-mterferometry and direct imaging to find companions 
in the separation range 0.13" to 6". After correction for chance alignment with background 
stars, we find a multiplicity (number of binaries or multiples divided by number of systems) of 
(14.0 ± 4.3) % and a companion star frequency (number of companions divided by number of 
systems) of (14.7 ± 5.1)%. Compared to solar-type main-sequence stars, the companion star 
frequency is lower by a factor of 0.61 ± 0.27. This is remarkably different from the high multiplic- 
ity found in the Taurus- Auriga star-forming region and for T Tauri stars in Chamaeleon known 
before ROSAT. We find only a few binaries with projected separations of more than 70 AU, also 
in contrast to the results for stars known before ROSAT. This indicates that the X-ray selected 
stars belong to a different population than the stars known before ROSAT, a hypothesis further 
supported by their Hipparcos distances and proper motions. 
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Recent studies of multiplicity in several star- 
forming regions have shown that multiplicity is 
- besides the initial mass function - one of the 
key characteristics of the star-formation process. 
Surveys of T Tauri stars (TTS) in the T associa- 
tion Taurus- Auriga (Leinert et al. 1993; Ghez et al. 
1993; Kohler & Leinert 1998) found an overabun- 
dance of binaries compared to solar-type main- 
sequence stars (Duquennoy & Mayor 1991) by a 
factor of two. On the other hand, surveys of 
low-mass stars in the Orion Trapezium Cluster 
(Prosser et al. 1994; Padgett et al. 1997; Petr et 
al. 1998) found a multiplicity that is compara- 
ble to or even lower than among main-sequence 
stars. Scally et al. (1999) used proper motion 
data of stars in the Trapezium cluster to search 
for binaries with separations in the range 1000 



1 Based on observations obtained at the European 
Southern Observatory, La Silla, proposal number 56.E-0197 



to 5000 AU. Their statistical analysis shows that 
the cluster contains no binaries at all in this sep- 
aration range. Studies of other clusters (IC348: 
Duchene et al. 1999; Pleiades: Bouvier et al. 1997) 
also found a binary frequency comparable to main- 
sequence stars. An explanation for these results 
is the disruption of binaries in stellar encoun- 
ters, which happen predominantly in dense clus- 
ters (Kroupa 1995b). However, it is also possi- 
ble that the multiplicity depends on environmen- 
tal parameters like magnetic fields or the cloud 
temperature (e.g. Durisen & Sterzik 1994). 

Kroupa (1995a) proposed to use the binary fre- 
quency and their distribution of periods to deter- 
mine the dominant mode of star formation, i.e. 
if most stars form in dense regions like the Orion 
Cluster, or in loose associations like Taurus- Auriga 
("inverse population synthesis"). However, al- 
though Taurus- Auriga is the best studied T asso- 
ciation, it is by no means clear that it is a typical 
example with respect to the formation of binaries. 
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Fig. 1. — Spatial distribution of the stars in our sample and the sample of Ghez et al. (1997). The contours 
mark the emission found in the IRAS 100 /im survey. Cha I and II are the cloud complexes where most TTS 
known before ROSAT are located. The large dashed circle encloses the four stars in the 77 Chamaeleontis 
cluster (Mamajek et al. 1999). 



Some authors surveyed other star-forming regions 
(Rcipurth & Zinnecker 1993; Brandner et al. 1996; 
Ghez et al. 1997), but the number of stars ob- 
served or the spatial resolution was very limited, 
which causes large statistical errors. Furthermore, 
it turned out that even within one star-forming re- 
gion (Scorpius-Centaurus), the period distribution 
of binaries in two subregions can be significantly 
different (Brandner & Kohler 1998; Kohlcr et al. 
2000). We can expect to find even larger differ- 
ences between different star-forming regions. 

In the light of these results, it is desirable to 
conduct a survey of a large number of stars in 
other T associations, in order to find out if Taurus- 
Auriga is a typical example or the exception to 
the rule. Thanks to the ROSAT All-Sky Survey, 
sufficiently large samples of TTS are available for 



the T associations Chamaeleon and Lupus. We 
carried out multiplicity surveys for both of them. 
This paper describes the survey in Chamaeleon, 
Lupus will be presented elsewhere. 

Ghez et al. (1997) already searched for bina- 
ries among TTS in Chamaeleon. However, only 
the observations of 17 stars were sensitive enough 
to give meaningful results. Furthermore, they ob- 
served only stars known before ROSAT, i.e. stars 
located on or near the dark clouds and mainly 
Classical T Tauri stars (CTTS). The results of the 
ROSAT All-Sky survey have shown that prc-main- 
scquence stars can be found in much larger ar- 
eas around star-forming regions, and that the ma- 
jority of prc-main-scquence stars are Weak- line T 
Tauri stars (WTTS; Walter et al. 1988; Neuhauser 
et al. 1995). 
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Our survey encompasses a list of WTTS found 
with the help of ROSAT. Section 2 describes the 
object list and section 3 our observations and data 
reduction. The results are presented and discussed 
in section 3 and 4, while section 5 contains the 
summary and conclusions. 

2. The Sample 

Our survey is based on the results of Alcala 
ct al. (1995) and Covino ct al. (1997). Alcala et 
al. did follow-up observations of optical counter- 
parts of X-ray sources in the Chamaeleon region 
detected in the ROSAT All-Sky Survey. They 
identified 82 young low-mass stars based on the 
presence of the Li A6707 absorption line and the 
spectral type. A complete list with coordinates, 
photometric data, and finding charts can be found 
in Alcala et al. (1995). Fig. 1 shows the spatial dis- 
tribution of these stars and the stars observed by 
Ghcz ct al. (1997). 

The four stars around 8 h 42 m , —79° are proba- 
bly members of the r\ Chamaeleontis cluster (Ma- 
majek et al. 1999). Our results for these stars are 
listed in Table 1 and 2, but we do not use them 
for the discussion of the binary properties in the 
Chamaeleon star-forming region. The binaries in 
j] Cha will be discussed in Kohler (2001). Two 
sources (RXJ 1039.5-7538N+S) are separated by 
only 5", therefore we count them as one binary. 
This leaves us with a sample of 82 — 4 — 1 = 77 
systems. 

It has been questioned recently (Briceno et 
al. 1997) whether stars found with low-resolution 
spectroscopy are indeed pre-main-sequence (PMS) 
objects or if they are somewhat older stars which 
have already reached the zero-age main-sequence 
(ZAMS). In order to unambiguously distinguish 
PMS and ZAMS objects, Covino et al. (1997) car- 
ried out observations with high spectral resolution 
of some 70 stars of the sample. They compared the 
equivalent width of the Li A6707 line to that of 
Pleiades stars of the same spectral type and clas- 
sified the stars as PMS or ZAMS. They found that 
more than 50 % of the stars could be confirmed to 
be PMS objects. 

3. Observations and Data Reduction 

The speckle observations were carried out at the 
ESO 3.5 m New Technology Telescope (NTT) on 



RXJ 0933.2-7433 




Fig. 2. — Examples of our speckle data. The im- 
ages in the first row show modulus (left) and phase 
(right) of the complex visibility of RXJ 0933.2- 
7433, a close, but well-resolved binary. The central 
maximum and the first maximum on both sides are 
clearly detected. The second row shows modulus 
and phase of RXJ 1014.2-7636, a potential binary 
closer than the diffraction limit of the telescope. 
In this case, we see only the central maximum and 
the first minimum on each side. 

La Silla, Chile, from 29. March to 2. April 1996. 
We used the SHARP I camera (System for High 
Angular Resolution Pictures) of the Max-Planck- 
Institut for Extraterrestrial Physics (Hofmann et 
al. 1992). All observations were done in the K- 
band at 2.2 /im. 

Although speckle interferometry can be consid- 
ered by now a standard technique (Leinert 1992), 
no program for speckle data reduction was pub- 
licly available at the time this survey was started. 
Therefore we used our speckle program 2 , which 
was already used for the surveys in Taurus- Auriga 
(Kohler & Leinert 1998) and Scorpius-Centaurus 

2 Now available on the software web page of the Center for 
Adaptive Optics at 

http://www.ucolick.org/~cfao/distributedsw/ds.html and 
http://babcock.ucsd.edu/cfao_ucsd/software.html 
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(Kohler et al. 2000). In this program, the modu- 
lus of the complex visibility (i.e. the Fourier trans- 
form of the object brightness distribution) is de- 
termined from power spectrum analysis, the phase 
is computed using the Knox-Thompson algorithm 
(Knox & Thompson 1974), and from the bispec- 
trum (Lohmann et al. 1983). Figure 2 shows ex- 
amples of reconstructed complex visibilities. For a 
more detailed description see Kohler et al. (2000). 

If the object appears unresolved, we compute 
the maximum brightness ratio of a companion that 
could be hidden in the noise of the data. The 
principle is to determine how far the data devi- 
ate from the nominal result for a point source 
(modulus = 1, phase = 0). We then compute the 
maximum brightness ratio of a companion that 
would be compatible with this amount of devia- 
tion. This is repeated for different position angles, 
and the maximum is used as upper limit for the 
brightness ratio of an undetected companion. See 
Leinert et al. (1997) for a more detailed descrip- 
tion of this procedure. 

If the object is a binary or triple, we com- 
pute a multidimensional least-squares fit using 
the amoeba algorithm (Press et al. 1994). Our 
program tries to minimize the difference between 
modulus and phase computed from a model bi- 
nary or triple and the observational data by vary- 
ing the separation, position angle, and brightness 
ratio of the model. This is necessary because the 
reconstructed images are a complex function of 
the 2-dimensional separation vector and flux ra- 
tio that cannot be solved to compute the binary 
parameters directly from the data. Fits to dif- 
ferent subsets of the data yield an estimate for 
the standard deviation of the binary parameters. 
We then subtract the contribution of the compan- 
ion^) from the images and apply the procedure 
described in the previous paragraph to find limits 
for the brightness of an undetected companion. 

In order to find binaries that are separated by 
more than 3", we obtained additional infrared 
images with the ESO/MPIA 2.2 m telescope on 
La Silla in February and March 1996 using the 
IRAC2b camera. For data reduction, we used the 
Daophot package within IRAF 3 . 

3 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under coop- 



4. Results 

4.1. Uncorrected Data 

All the binary and multiple stars we find in our 
sample are listed in Table 1, while Table 2 gives 
limits for the brightness of undetected compan- 
ions. Figure 3 shows these results as a plot of flux 
ratio and magnitude difference vs. binary separa- 
tion. In total, among the 81 systems of the sample 
we find 18 binaries and one triple star. 

The modulus of the complex visibility of a bi- 
nary is a cosine-shaped function. If the separation 
of the binary is equal to the diffraction limit (0.13" 
for a 3.5m telescope at K), exactly one period of 
the modulus of the visibility fits within the radius 
where the optical transfer function of the telescope 
is not zero. Under good circumstances, it is pos- 
sible to discover binaries with even smaller sepa- 
rations, down to about half the diffraction limit 
(Table 1 shows that we actually do find some, an 
example is presented in Figure 2). However, in 
these cases we can detect only the first minimum, 
but not the second maximum of the modulus of 
the visibility. Therefore, we cannot distinguish 
with certainty a close binary star from an elon- 
gated structure. Even more important, we can- 
not be sure that we find all companions at sepa- 
rations less than the diffraction limit. For these 
reasons, we limit ourselves to companions in the 
separation range between 0.13" and 6". The up- 
per limit was chosen so that contamination with 
background stars has little effect (see section 4.3 
for a detailed discussion of this problem). 

We found four binaries with separations smaller 
than 0.13". These binaries are treated as unre- 
solved stars in our survey. This leaves 14 binaries 
and one triple, i. e. 16 companion stars, including 
one binary in the rj Cha cluster. For ease of com- 
parison we mention that these uncorrected data 
correspond to a fractional multiplicity (number of 
multiples divided by total number of systems) of 
0.19 ± 0.05 or to a number of 0.20 ± 0.05 com- 
panions per primary. These numbers change to 
0.18 ± 0.05 and 0.19 ± 0.05 if we exclude the four 
stars in the rj Chamaeleontis cluster. 



erative agreement with the National Science Foundation. 
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Fig. 3. — The results of our multiplicity survey in a plot of flux ratio or magnitude difference vs. binary star 
separation. The data points mark the detected companion stars, the thick line shows the average and the thin 
line the worst sensitivity for undetected companions. The dashed vertical line at 0.13" shows the diffraction 
limit for a 3.5 m telescope at K. This is the limit for unambiguous identification of binary stars. The dashed 
horizontal line shows the completeness limit in flux ratio for the whole survey. The two observations of 
RXJ 1039.5-7538 are represented by two data points connected by a line. 



4.2. Completeness 

Figure 3 shows not only the stars where we find 
companions, but also the sensitivity of our survey, 
i.e. the maximum brightness ratio of a possible 
undetected companion as a function of the sepa- 
ration. This sensitivity depends on factors like the 
atmospheric conditions at the time of the observa- 
tions and the brightness of the target star. Fig. 3 
and Table 2 show that in most cases we are able 
to detect companions with magnitude differences 
to the primary up to (at least) 2.5 mag . This cor- 
responds to a brightness ratio of 1:10. 

Our observations of four stars are not suffi- 
ciently sensitive to exclude faint companions at 



small separations. However, based on the number 
of companions actually found, we expect about 
0.3 additional companions above a brightness ra- 
tio of 1:10 at separations > 0.13". Therefore, we 
are confident we have found all companions with 
separations between 0.13" and 6" and with a mag- 
nitude difference of less than 2.5 mag . 
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Table 1 
Binary and triple stars. 



No. Designation 


Evol. 


Date of Separation 


Position Brightness 




Status a 


Observation 


Angle [°] Ratio at K 



1 


RXJ 


0837. 


0-7856 b 


PMS 


29. Mar. 


96 


0. 


135±0, 


.003 


15. 


.9±1. 


,8 


0. 


,805±0. 


,024 


11 


RXJ 


0915. 


.5-7609 


PMS 


30. Mar. 


96 


0. 


lllzbO. 


,007 


292 


,5±4. 


.3 


0. 


,255±0. 


,02 


13 


RXJ 


0919 


.4-7738N 


? 


30. Mar. 


96 


0. 


109±0. 


,003 


173 


,9±1, 


,2 


0. 


658±0. 


,021 


16 


RXJ 


0933 


.2-7433 


ZAMS 


31. Mar. 


96 


0. 


,222±0. 


,003 


231 


,4±0, 


,4 


0. 


,743±0. 


,034 


17 


RXJ 


0935 


,0-7804 


PMS 


31. Mar. 


96 


0, 


,360±0. 


,003 


353 


,9±0. 


,2 





,901±0. 


,022 


22 


RXJ 


0952. 


.7-7933 


? 


31. Mar. 


96 


0, 


,267±0. 


,003 


314 


.6±0. 


.5 





,259±0. 


,018 


26 


RXJ 


1009 


,6-8105 


? 


28. Feb. 


96 


5, 


,413±0 


,017 


136 


,9±0. 


,3 





,039±0. 


,001 


27 


RXJ 


1014. 


,2-7636 


? 


1. Apr. 


96 


0. 


,091±0, 


,007 


259 


,6±6. 


.5 


0. 


,593±0, 


,062 


31/32 


RXJ 


1039 


.5-7538N+S 


? 


26. Feb. 


96 


4, 


,962±0. 


.022 


347 


,4±0. 


.2 


0. 


,287±0, 


003 












30. Mar. 


96 


5, 


,002±0. 


.003 


347 


,5±0, 


.1 


0. 


,264±0, 


,004 


36 


RXJ 


1108. 


,8-7519b 


PMS 


30. Mar. 


96 


0, 


,150±0. 


.003 


7. 


,8±1, 


.7 


0, 


,716±0. 


.019 


43 


RXJ 


1125, 


.8-8456 


? 


26. Feb. 


96 


3, 


,545±0. 


.022 


208 


,0±0. 


.2 


0, 


140±0. 


.003 


44 


RXJ 


1129. 


.2-7546 


PMS 


29. Mar. 


96 


0, 


534±0. 


.003 


28. 


,7±0. 


.2 


0, 


189±0. 


.008 


48 


RXJ 


1150, 


.9-7411 


? 


1. Apr. 


96 


0, 


,875±0. 


.003 


106 


,0±0. 


,1 


0, 


,257±0. 


.012 


50 


RXJ 


1158. 


5-7754a 


PMS 


29. Mar. 


96 


0, 


,073±0. 


,006 


160 


,8±2. 


.8 


0, 


.506±0 


.089 


57 


RXJ 


1203 


.7-8129 


? 


31. Mar. 


96 


0. 


,396±0 


.003 


238 


,3±0. 


.1 





,775±0. 


,008 


60 


RXJ 


1207. 


.9-7555 


? 


1. Apr. 


96 


0. 


,588±0. 


.003 


10. 


,2±0. 


.2 





,924±0, 


,02 


65 


RXJ 


1220 


.4-7407 


PMS 


1. Apr. 


96 


0. 


,296±0 


.003 


348 


,4±1. 


.2 





,258±0. 


,021 


73 


RXJ 


1243 


.1-7458 A-B 


PMS 


30. Mar. 


96 


0, 


,280±0. 


.004 


84. 


,7±3. 


.6 





,386±0. 


.025 








AB-C 




30. Mar. 


96 


2, 


,513±0 


.003 


258 


,7±0. 


.1 





,244±0. 


,007 


75 


RXJ 


1301. 


.0-7654 


PMS 


30. Mar. 


96 


1, 


,444±0, 


.005 


3 


,2±0, 


.2 


0. 


,155±0, 


.01 



a Covino ct al. (1997) 

b member of the r\ Chamaeleontis Cluster (Mamajek et al. 1999) 
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Table 2 

Limits for the brightness of undetected companions. 



No. Designation Evol. Date of Maximum Flux Ratio Minimal Atok [mag] 









Status a 


Observation 


at 0.13" 


at 0.5" 


at 0.13" 


at 0.5' 


1 


RXJ 


0837.0-7856 b 


PMS 


29. 


Mar. 


96 


0.04 


0.01 


3.49 


5.00 


2 


RXJ 


0842.4-8345 


? 


30. 


Mar. 


96 


0.05 


0.02 


3.25 


4.25 


3 


RXJ 


0842.9-7904 b 


PMS 


29. 


Mar. 


96 


0.05 


0.01 


3.25 


5.00 


4 


RXJ 


0844.5-7846 b 


PMS 


29. 


Mar. 


96 


0.05 


0.02 


3.25 


4.25 


5 


RXJ 


0848.0-7854 b 


PMS 


29. 


Mar. 


96 


0.05 


0.02 


3.25 


4.25 


6 


RXJ 


0849.2-7735 


? 


29. 


Mar. 


96 


0.04 


0.01 


3.49 


5.00 


7 


RXJ 


0850.1-7554 


PMS 


29. 


Mar. 


96 


0.06 


0.01 


3.05 


5.00 


8 


RXJ 


0853.1-8244 


? 


30. 


Mar. 


96 


0.05 


0.02 


3.25 


4.25 


9 


RXJ 


0901.0-7715 


? 


30. 


Mar. 


96 


0.2 


0.07 


1.75 


2.89 


10 


RXJ 


0902.9-7759 


PMS 


30. 


Mar. 


96 


0.05 


0.02 


3.25 


4.25 


11 


RXJ 


0915.5-7609 


PMS 


30. 


Mar. 


96 


0.05 


0.01 


3.25 


5.00 


12 


RXJ 


0917.2-7744 


? 


30. 


Mar. 


96 


0.07 


0.01 


2.89 


5.00 


13 


RXJ 


0919.4-7738N 


? 


30. 


Mar. 


96 


0.04 


0.01 


3.49 


5.00 


14 


RXJ 


0919.4-7738S 


? 


30. 


Mar. 


96 


0.17 


0.07 


1.92 


2.89 


15 


RXJ 


0928.5-7815 


? 


30. 


Mar. 


96 


0.08 


0.02 


2.74 


4.25 


16 


RXJ 


0933.2-7433 


ZAMS 


31. 


Mar. 


96 


0.06 


0.01 


3.05 


5.00 


17 


RXJ 


0935.0-7804 


PMS 


31. 


Mar. 


96 


0.04 


0.01 


3.49 


5.00 


18 


RXJ 


0936.3-7820 


? 


31. 


Mar. 


96 


0.06 


0.02 


3.05 


4.25 


19 


RXJ 


0942.7-7726 


PMS 


31. 


Mar. 


96 


0.03 


0.01 


3.81 


5.00 


20 


RXJ 


0946.9-8011 


? 


31. 


Mar. 


96 


0.05 


0.03 


3.25 


3.81 


21 


RXJ 


0951.9-7901 


PMS 


31. 


Mar. 


96 


0.04 


0.02 


3.49 


4.25 


22 


RXJ 


0952.7-7933 


? 


31. 


Mar. 


96 


0.05 


0.01 


3.25 


5.00 


23 


RXJ 


1001.1-7913 


PMS 


31. 


Mar. 


96 


0.05 


0.01 


3.25 


5.00 


24 


RXJ 


1005.3-7749 


PMS 


1 


. Apr. 


96 


0.06 


0.01 


3.05 


5.00 


25 


RXJ 


1007.7-8504 


? 


1 


Apr. 


96 


0.05 


0.01 


3.25 


5.00 


26 


RXJ 


1009.6-8105 


? 


1 


.Apr. 


96 


0.06 


0.02 


3.05 


4.25 


27 


RXJ 


1014.2-7636 


? 


1 


.Apr. 


96 


0.05 


0.01 


3.25 


5.00 


28 


RXJ 


1014.4-8138 


? 


1 


. Apr. 


96 


0.04 


0.02 


3.49 


4.25 


29 


RXJ 


1017.9-7431 


ZAMS 


1 


.Apr. 


96 


0.05 


0.02 


3.25 


4.25 


30 


RXJ 


1035.8-7859 


ZAMS 


30. 


Mar. 


96 


0.06 


0.03 


3.05 


3.81 


31/32 


RXJ 


1039.5-7538N+S 


? 


30. 


Mar. 


96 


0.12 


0.05 


2.30 


3.25 


33 


RXJ 


1044.6-7849 


ZAMS 


30. 


Mar. 


96 


0.10 


0.02 


2.50 


4.25 


34 


RXJ 


1048.9-7655 


? 


30. 


Mar. 


96 


0.06 


0.03 


3.05 


3.81 


35 


RXJ 


1108.8-7519a 


PMS 


30. 


Mar. 


96 


0.13 


0.04 


2.22 


3.49 


36 


RXJ 


1108.8-7519b 


PMS 


30. 


Mar. 


96 


0.05 


0.01 


3.25 


5.00 


37 


RXJ 


1109.4-7627 


PMS 


30. 


Mar. 


96 


0.09 


0.02 


2.61 


4.25 


38 


RXJ 


1111.7-7620 


PMS 


30. 


Mar. 


96 


0.06 


0.03 


3.05 


3.81 


39 


RXJ 


1112.7-7637 


PMS 


29. 


Mar. 


96 


0.07 


0.05 


2.89 


3.25 


40 


RXJ 


1117.0-8028 


PMS 


30. 


Mar. 


96 


0.05 


0.03 


3.25 


3.81 


41 


RXJ 


1120.3-7828 


? 


29. 


Mar. 


96 


0.04 


0.02 


3.49 


4.25 


42 


RXJ 


1123.2-7924 


PMS 


30. 


Mar. 


96 


0.07 


0.05 


2.89 


3.25 
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Table 2 — Continued 



No. Designation Evol. Date of Maximum Flux Ratio Minimal Atok [mag] 











Status a 


Observation 


at 1 3" 


at 0.5" 


at 1 3" 


at 0.5' 


43 


RXJ 


1125 


.8-8456 


? 


OU. 


A A Q T* 




0.05 


0.02 


3.25 


4.25 


44 


RXJ 


1129 


.2-7546 


PMS 


90 
zy . 


IV /To r 




0.05 


0.01 


3.25 


5.00 


45 


RXJ 


1140 


.3-8321 


? 


OU. 


A T q r 


yu 


0.03 


0.03 


3.81 


3.81 


46 


RXJ 


1149 


.8-7850 


PMS 


1 

1 




yo 


0.04 


0.03 


3.49 


3.81 


47 


RXJ 


1150 


.4-7704 


PMS 


90 

zy . 


A/To i- 

rviar. 


yo 


0.10 


0.02 


2.50 


4.25 


48 


RXJ 


1150 


.9-7411 


? 


i 
i 


. Apr. 


yo 


0.08 


0.01 


2.74 


5.00 


49 


RXJ 


1157 


.2-7921 


PMS 


90 

zy . 


A/To 1- 

iviar. 


yo 


0.10 


0.03 


2.50 


3.81 


50 


RXJ 


1158 


.5-7754a 


PMS 


90 

zy . 


A/To-r 

iviar. 


yo 


0.03 


0.01 


3.81 


5.00 


51 


RXJ 


1158 


,5-7754b 


PMS 


90 

zy . 


A/Tq r 


yo 


0.05 


0.03 


3.25 


3.81 


52 


RXJ 


1158 


.5-7913 


PMS 


90 
zy . 


A/To r 


yo 


0.06 


0.01 


3.05 


5.00 


53 


RXJ 


1159 


.7-7601 


PMS 


1 . 


A/Tot- 

IVltil . 


yo 


0.02 


0.01 


4.25 


5.00 


54 


RXJ 


1201 


.7-7859 


PMS 


01. 


A/To v 

IVltil . 


yo 


0.06 


0.01 


3.05 


5.00 


55 


RXJ 


1202 


.1-7853 


PMS 


01. 


A/To-r 


yo 


0.04 


0.01 


3.49 


5.00 


56 


RXJ 


1202 


.8-7718 


? 


01. 


A T Q T* 

IVldil . 


yo 


0.03 


0.01 


3.81 


5.00 


57 


RXJ 


1203 


.7-8129 


? 


(31. 


A T o t* 

iviar . 


yo 


0.03 


0.01 


3.81 


5.00 


58 


RXJ 


1204 


.6-7731 


PMS 


Ol. 


A/To 1- 
IVldX. 


yo 


0.04 


0.01 


3.49 


5.00 


59 


RXJ 


1207 


.7-7953 


? 


1 
1 


A nr 

. Apr. 


GR 

yo 


0.07 


0.03 


2.89 


3.81 


60 


RXJ 


1207 


.9-7555 


? 


1 
1 


A nr 

. Apr. 


Gfi 

yo 


0.04 


0.01 


3.49 


5.00 


61 


RXJ 


1209 


.8-7344 


? 


1 
1 


. Apr. 


Gfi 

yo 


0.04 


0.02 


3.49 


4.25 


62 


RXJ 


1216 


.8-7753 


PMS 


1 

1 


. Apr. 


Gfi 

yo 


0.10 


0.04 


2.50 


3.49 


63 


RXJ 


1217 


.4-8035 


? 


1 

1 


. Apr. 


yo 


0.07 


0.01 


2.89 


5.00 


64 


RXJ 


1219 


.7-7403 


PMS 


1 

1 


. Apr. 


yo 


0.04 


0.02 


3.49 


4.25 


65 


RXJ 


1220 


.4-7407 


PMS 


1 

1. 


. Apr. 


yo 


0.05 


0.01 


3.25 


5.00 


66 


RXJ 


1220 


.6-7539 


? 


1 

1 


. Apr. 


yo 


0.03 


0.01 


3.81 


5.00 


67 


RXJ 


1223 


.5-7740 


? 


1 

1 


. Apr. 


yo 


0.04 


0.02 


3.49 


4.25 


68 


RXJ 


1224. 


.8-7503 


? 


1 
1 


. Apr. 


Gfi 

yo 


0.04 


0.02 


3.49 


4.25 


69 


RXJ 


1225. 


.3-7857 


? 


1 

1 


. A_pr. 


yu 


0.08 


0.01 


2.74 


5.00 


70 


RXJ 


1231 


.9-7848 


? 


1 


. Apr. 


96 


0.06 


0.03 


3.05 


3.81 


71 


RXJ 


1233 


.5-7523 


? 


30. 


Mar. 


96 


0.04 


0.02 


3.49 


4.25 


72 


RXJ 


1239 


.4-7502 


PMS 


29. 


Mar. 


96 


0.05 


0.01 


3.25 


5.00 


73 


RXJ 


1243 


.1-7458 


PMS 


30. 


Mar. 


96 


0.1 


0.02 


2.50 


4.25 


74 


RXJ 


1243 


.6-7834 


? 


30. 


Mar. 


96 


0.10 


0.03 


2.50 


3.81 


75 


RXJ 


1301 


.0-7654 


PMS 


30. 


Mar. 


96 


0.09 


0.02 


2.61 


4.25 


76 


RXJ 


1303 


.3-7706 


? 


30. 


Mar. 


96 


0.08 


0.04 


2.74 


3.49 


77 


RXJ 


1307 


.3-7602 


? 


30. 


Mar. 


96 


0.05 


0.04 


3.25 


3.49 


78 


RXJ 


1320 


.0-7406 


? 


31. 


Mar. 


96 


0.11 


0.05 


2.40 


3.25 


79 


RXJ 


1325. 


.7-7955 


? 


31. 


Mar. 


96 


0.06 


0.01 


3.05 


5.00 


80 


RXJ 


1346 


.4-7409 


? 


31. 


Mar. 


96 


0.12 


0.04 


2.30 


3.49 


81 


RXJ 


1349 


.2-7549 


? 


1 


.Apr. 


96 


0.03 


0.01 


3.81 


5.00 


82 


RXJ 


1415. 


.0-7822 


PMS 


1 


.Apr. 


96 


0.06 


0.01 


3.05 


5.00 



a Covino ct al. (1997) 

b member of the r] Chamaeleontis cluster (Mamajek et al. 1999) 
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4.3. Confusion with Background Stars 



Since we observed our targets only once, we 
have no possibility to detect orbital motion. 
Therefore, we cannot say if a given binary is in- 
deed a physically bound pair or only a chance 
alignment with a field star. To estimate the num- 
ber of chance projections, we count the field stars 
in 76 of the infrared images obtained at the ESO / 
MPIA 2.2 m telescope. We exclude a circular area 
with a radius of 15" around the T Tauri star in 
each image to avoid counting physically bound 
companions. This leaves 4400arcsec 2 per field, 
and a total area of 93 arcmin 2 . 

The distribution of the number of field stars we 
obtain in this way is shown in Fig. 4. Ovcrplotted 
is a Poisson distribution with a mean of 1.89, the 
average number of stars found per field. The Pois- 
son distribution gives a reasonable fit to the data, 
indicating that the number of field stars can be 
described by Poisson statistics. The average num- 
ber of stars per field yields a field star density of 
(4.3 ± 0.5) • 10~ 4 stars per arcsec 2 . So the proba- 
bility to find a field star with a projected distance 
of less than 6" to one T Tauri star is 

(4.3 ± 0.5) • 10~ 4 • 7T • 6 2 = (4.9 ± 0.6)%. 

The expected number of chance-projected field 
stars in our sample of 77 targets is 

(4.9 ±0.6)% -77 = 3.7 ±0.4. 

The number of physically bound companions is 
therefore 15 — 3.7 w 11 companion stars. To obtain 
the number of bound binaries and triples, we must 
take into account that the number of "false" bina- 
ries depends on the number of "true" single stars, 
and the number of "false" triples on the number 
of "true" binaries. The result of this trivial, but 
somewhat tedious calculation is 10.3 ± 0.4 bound 
binaries and 0.5 ± 0.6 bound triples. This means 
we cannot tell if the triple system RXJ 1243.1-7458 
consists of three bound TTS or a binary TTS and 
a field star. However, the probability for a chance 
projection is proportional to the separation. Since 
the pair AB-C of the triple is one of the four widest 
pairs in our sample, component C is probably a 
field star. 

These corrected numbers correspond to a frac- 
tional multiplicity of (10.3 + 0.5)/77 = 0.140 ± 
0.043 and to 11.3/77 = 0.147 ± 0.051 companions 



50 o • 



O 

d o 



I 1 1 , 1 1 1 , 1 1 1 , 1 1 . 1 1 . 

3 4 6 8 

Number of Field Stars per Image 

Fig. 4. — Distribution of the number of field stars 
in 76 images. The histogram shows the number 
of background stars we count in our images; the 
dots denote a Poisson distribution with the same 
average star density 

per primary. The stars in the r\ Chamaeleontis 
cluster have been excluded to obtain this result. 

4.4. Bias Induced through X-Ray Selec- 
tion 

Brandner ct al. (1996) pointed out that the flux 
limit of X-ray selected samples induces a detec- 
tion bias in favour of binary and multiple systems. 
There is a small number of systems with a com- 
bined X-ray flux above the detection threshold, al- 
though the fluxes from individual components are 
below the cut-off. These systems cause an overes- 
timate of the multiplicity. 

To check if our results are affected by this bias, 
we look at the X-ray counts that led to the de- 
tection of the stars (Alcala ct al. 1995) in order 
to search for binaries that would have been de- 
tected only because of the bias. These systems 
would have a number of counts below twice the 
detection threshold, since otherwise at least one 
component would be bright enough to be detected 
without the additional X-ray flux from its com- 
panion. This is a very conservative limit since not 
all binaries consist of stars with equal X-ray fluxes. 

The distributions of X-ray counts of unresolved, 
binary and triple stars are plotted in Fig. 5. The 
lowest number of counts measured is 12.5, while 
the usually adopted value for the detection limit 
of the RASS is 8. In Fig. 5, we marked 10 counts 
in the panel for unresolved stars, 20 counts in the 
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Fig. 6. — Simulated distribution of projected sep- 
arations for a sample of 10 million main-sequence 
binaries with a system mass of 1M Q . The his- 
234 togram shows the simulated data, the line is a 

log-normal distribution fitted to the histogram. 

log 10 (Counts) 



Fig. 5. — Distribution of the number of X-ray 
counts of unresolved, binary, and triple systems. 
The vertical lines mark the limits chosen by us to 
obtain an unbiased sample: Mimit = 10 counts for 
unresolved stars, 20 counts for binaries, and 30 
counts for triple stars 

binary panel, and 30 counts in the panel for the 
triple. The figure shows that none of our multi- 
ple systems falls even close to the detection limit. 
Only two unresolved stars could be affected if they 
would be binaries with companions undetectable 
by our survey. However, since the binary fraction 
of our sample is below 20 %, we have no reason to 
assume that this is the case. 

We conclude that we do not have to apply any 
corrections for a bias induced by the X-ray selec- 
tion of our sample. 

5. Discussion 

5.1. The Distribution of Separations of 
Main- Sequence Stars 

In order to compare the young binaries to 
main-sequence stars, we have to face the prob- 
lem that we observed the stars only once, so we 
know separation and position angle at one date. 
(RXJ 1039.5-7538 was observed twice, but the 
timespan of one month is not long enough to de- 
tect orbital motion.) Duqucnnoy & Mayor (1991) 



give the period distribution of main-sequence bi- 
naries, which we cannot derive from our data. In 
preceding papers (e.g. Leinert et al. 1993; Kohlcr 
& Leinert 1998; Kohler et al. 2000), we used statis- 
tical arguments to convert the separation distribu- 
tion of binaries to a period distribution. This relies 
on assumptions on the distributions of orbital el- 
ements like eccentricity and inclination. Further- 
more, given the small number of multiple stars in 
our sample, it is not clear how reliable the result 
of this statistical argumentation is. 

In this work, we go the other way and use the 
well-known distribution of main-sequence binaries 
to compute their distribution of projected separa- 
tions. To do this, we simulated samples of typi- 
cally 10 million binaries with orbital elements dis- 
tributed according to Duqucnnoy & Mayor (1991), 
i.e. the periods have a log-normal distribution, 
the distribution of eccentricities is f(e) = 2e, 
the inclinations are distributed isotropically, while 
all the other parameters follow uniform distribu- 
tions. In most of the simulations, we kept the 
system mass fixed, but we performed simulations 
with different system masses between 0.5 M© and 
2 M Q . We also simulated samples where the sys- 
tem mass is not fixed, but distributed uniformly 
in the ranges 0.5 .. . 1.5 M Q , 0.8 .. . 1.5 M Q , and 
0.5 ... 2 M©. The mass ratio of the binaries does 
not enter the computation, so we do not have to 
worry about its distribution. 
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Figure 6 shows a typical result of these simu- 
lations. The projected separations follow a log- 
normal distribution of the form 

/(aO = C-exp(-fc?^), 

where x is the logarithm of the projected sepa- 
ration in AU. To obtain the constants x and a, 
we fitted a Gaussian function to the simulated 
data. The result for x depends slightly on the 
system mass and varies between 1.37 and 1.51, 
while all the simulations resulted in the same a 
of 1.55. In the following, we use x = 1.44, the re- 
sult of the simulation with a constant system mass 
of 1 M 4 . To get the constant C, we use the num- 
ber of companions given by Duqucnnoy & Mayor 
(1991). They found 101 companions in their sam- 
ple of 164 main-sequence stars, so we adjust C 
such that the integral over f(x) is 101/164. 

According to the conversion from projected sep- 
aration to period used in previous papers, a sep- 
aration of 10 144 AU corresponds to a period of 
10 days, and o\ ogs — 1.55 translates to <J\ gP = 
2.3. These are the parameters of the period dis- 
tribution given by Duquennoy & Mayor (1991), 
which demonstrates that the old method was ap- 
propriate to compare a sample of young binaries to 
main-sequence stars in the statistical sense. How- 
ever, it gave the false impression that we are able 
to measure periods of binaries based on a single ob- 
servation. Therefore, we will use the new method 
in the following and compare the distributions of 
separations, not periods. 

5.2. The Distance to Chamaeleon 

The final step needed before we can compare 
our results to those of Duquennoy & Mayor (1991) 
is to convert the projected separations from arc- 
seconds to AU. To do so, we need the distance to 
the binary stars. Hipparcos parallaxes of three 
stars in Cha I give a mean distance of 171 ± 
20 pc (Wichmann et al. 1998), in agreement with 

4 Alcala et al. (1997) placed the ROSAT-selected stars in the 
Chamaeleon association in the Hcrtzsprung-Russcl diagram 
in order to determine their masses and ages. They found 
masses in the range 0.2 — 2.5 Mq with about half of the 
stars having masses below IMq. However, they assumed 
a distance of 150 pc. If the stars are actually closer (see 
Sect. 5.2), their real masses will be somewhat smaller, in 
the range 0.15 - 1.7 M 



the distance of about 140 - 150 pc that was ob- 
tained by various other methods (sec Schwartz 
(1991) for a review). Hipparcos measurements of 
ROSAT-discovered stars in Chamaeleon resulted 
in much smaller distances between 63 and 128 pc 
(Neuhauser & Brandner 1998). However, four of 
the seven stars observed by Hipparcos are bina- 
ries, which makes their parallax unreliable. The 
average distance of the other three stars is 65 pc 
if we use the errors to weight the individual dis- 
tances, and 74 pc if we do not take the errors into 
account. 

Frink et al. (1998) studied the kinematics of 
T Tauri stars in Chamaeleon, both stars known 
before and stars discovered by ROSAT. They 
found that their sample can be divided into at least 
two subgroups, where subgroup 1 contains mainly 
TTS known before ROSAT, while subgroup 2 con- 
sists mainly of ROSAT-discovered TTS. The mean 
proper motion of subgroup 2 is about twice as high 
as that of subgroup 1. If both subgroups have the 
same space velocity, this would imply that the 
mean distance of subgroup 1 is about twice the 
mean distance of subgroup 2, in agreement with 
the Hipparcos results. 

For a given star in our sample, it is usually not 
possible to decide if it belongs to the Chamae- 
leon association at a distance of about 140-170 pc, 
or to another group in the foreground at about 
70-80 pc. Fortunately, the difference in distance 
between 70 and 170 pc corresponds to a shift in 
log(s) of about 0.4, which is small compared to 
the width of the separation distribution of main- 
sequence stars. In the following, we will adopt 
a distance of 80 pc for the ROSAT-selected stars, 
and twice as much for the stars observed by Chez 
et al. (1997). 

5.3. Comparison to Main- Sequence Stars 

At a distance of 80 pc, the separation range 0.13 
- 6" corresponds to 10.4 - 480 AU. Our simula- 
tion of 1 M main-sequence binaries shows that 
(39.6±6.3) % of the companions are at a projected 
separation in this range. We adopted an error of 
v/to/101 = 6.3 %, because 40 is 39.6% of the 101 
binaries found by Duquennoy & Mayor (1991). Di- 
vided by the total number of 164 systems (single 
stars and multiples) in the main-sequence sample, 
this yields 24.1 ± 3.8 companions per 100 systems. 
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VI O 



Stars discovered by ROSAT 
/~\ Main — sequence stars 




12 3 

log(separation/AU) 

Fig. 7. — Companion star frequency as a func- 
tion of projected separation. The histogram shows 
the result of this work, the curve shows the distri- 
bution of separations for main-sequence binaries. 
This is the result of a simulation for system masses 
of 1 M©, as described in section 5.1. 



In the Chamaeleon sample of 77 systems, we 
find 11.3 ± 3.9 companions after correction for 
chance alignment with background stars. This 
corresponds to 14.7 ±5.1 companions per 100 sys- 
tems, which is 9.4 ± 6.3 less than among main- 
sequence stars. Thus, the companion star fre- 
quency in our sample is reduced by a factor of 
0.61 ±0.27 compared to main-sequence stars. This 
is surprising, since the Chamaeleon star-forming 
region is a T association similar to Taurus- Auriga, 
where we found a much higher binary frequency 
(Kohler & Leinert 1998). 

Fig. 7 shows the distribution of projected sep- 
arations of the stars in our sample and main- 
sequence stars. The companion star frequencies in 
the range 10 - 70 AU agree almost perfectly, but 
at larger separations, the number of companions 
to young stars drops to nearly zero and is signif- 
icantly lower than the number of companions to 
main-sequence stars. A Kolmogorov-Smirnov test 
shows that the two distributions are indeed differ- 
ent: the probability that both samples were drawn 
from the same distribution is only 6 • 10~ 6 . The 
separation distribution of young stars in Chamae- 
leon appears to be more similar to that of low-mass 
stars in the Orion Trapezium Cluster: a number 
of binaries with small separations that is compa- 



PMS stars 

Unknown evolutionary status 
Main— sequence stars 




12 3 

log(separation/AU) 

Fig. 8. — The distributions of projected separa- 
tions of PMS stars and stars of unknown evolu- 
tionary status. Only error bars for PMS stars are 
shown to enhance the clarity of the diagram. 

rable to main-sequence stars (Prosser et al. 1994; 
Padgett et al. 1997; Petr et al. 1998), and only a 
very small number of binaries at larger separations 
(Scally et al. 1999). 

5.4. Stars of different evolutionary status 

Covino et al. (1997) used high-resolution spec- 
tra to classify ROSAT-selected stars in Chamae- 
leon as PMS star, ZAMS star, or star of unknown 
evolutionary status. Their object list is only a sub- 
set of our sample; we treat stars not observed by 
them as stars of unknown evolutionary status. Ta- 
ble 3 shows the results of our multiplicity survey 
for the different classes. Due to the large errors 
caused by the small sample size, we are unable to 
find a significant difference between stars of differ- 
ent evolutionary status. 

Figure 8 shows the distribution of separations of 
PMS stars and stars of unknown evolutionary sta- 
tus (it would not be useful to plot the one ZAMS 
binary). The bias towards binaries with small 
separations appears to be even more pronounced 
among bona- fide PMS stars, while the distribu- 
tion of presumably older stars is more or less flat. 
We take this as an indication that the unusual 
distribution of separations we found for the whole 
sample is indeed a feature of the ROSAT-selected 
PMS stars in this region. 

A KS test gave a probability of 38 % that the 
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Table 3 
Companion star frequencies 





Main-sequence 


All Stars in 


PMS 


ZAMS 


Unknown 




Stars a 


Our Sample 






Evol. Status 


Total 


164 


77 


33 


4 


40 


Companions' 3 




11.3 ±3.9 


5.4 ±2.6 


0.8 ± 1.0 


5.0 ±2.6 


Companions per 100 systems 


24.1 ±3.8 


14.7 ±5.1 


16.3 ±8.0 


20.1 ±25.0 


12.6 ±6.6 



a Duquennoy & Mayor (1991) 

b projected separations in the range 10.4-480 AU, after subtraction of background stars 
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Fig. 9. — The distributions of projected separa- 
tions of stars observed by Ghez et al. (1997) and 
in this work. 

distributions of PMS stars and of stars of unknown 
evolutionary status were drawn from the same par- 
ent distribution. This result does not enable us to 
accept or reject the hypothesis that both distri- 
butions are identical. The sample of stars of un- 
known evolutionary status is probably a mixture 
of PMS and older stars. It shows a separation dis- 
tribution that is similar to main-sequence stars, 
albeit with a smaller total number of companions. 

5.5. Comparison to Ghez et al. (1997) 

Ghez et al. (1997) conducted a multiplicity 
speckle survey of TTS in southern star-forming 



regions known before ROSAT, including stars in 
the Chamaeleon star-forming region. The sensi- 
tivity of their observations for faint companions 
is rather non-uniform, but the observations of 17 
stars in Chamaeleon have a sensitivity compara- 
ble to our observations. We assumed a distance of 
160 pc for these stars, since the majority is located 
in the Cha I complex. 

Figure 9 shows the distribution of projected 
separations of the stars observed by Ghez et al. 
(1997) compared to the stars observed by us. The 
distributions are clearly different. A \ 2 test con- 
firms this, the probability that both samples were 
drawn from the same distribution is only 10 , 
a KS test gave a probability of 2.4%. This dif- 
ference cannot be explained by errors in the dis- 
tances. To bring the peaks of both distributions 
together, the distance to one of the groups would 
have to be changed by a factor of 10. 

This indicates that both samples represent two 
different populations of stars. This is not very sur- 
prising, since the distances of the two groups dif- 
fer by a factor of two (see section 5.2). However, 
observations of other star-forming regions have 
shown that regions with a higher stellar density 
contain fewer binaries (Petr et al. 1998; Duchcnc 
et al. 1999). In Chamaeleon, we find the opposite 
trend: the stars observed by Ghez et al. (1997) 
are located in the dense cores, yet their binary 
frequency is lower. 

There are two possible explanations for this 
result: either the dependence of the binary fre- 
quency on the environmental conditions is not a 
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simple function of stellar density, or the ROSAT- 
selected stars in Chamaeleon formed in a much 
denser environment. The proper motions of 31 of 
the ROSAT-sclccted stars are known (Frink et al. 
1998; Frink 1999). Their proper motion vectors 
are not oriented away from a single point, so they 
show no sign that the distribution was more con- 
centrated in the past. 

It is more likely that the stars formed near 
their present location in small "cloudlets" (Feigel- 
son 1996). Such small clouds have been discov- 
ered in a 13 CO (J = 1-0) survey by Mizuno et al. 
(1998). They found that about 70% of the TTS 
in the region surveyed are located within 4 pc of a 
small cloud. However, it is not clear why only a 
small number of binaries formed in the cloudlets 
towards Chamaeleon, while under similar condi- 
tions in Taurus- Auriga most stars formed in bi- 
nary systems. 

6. Summary and Conclusions 

We conducted a multiplicity survey of some 
80 X-ray selected young stars. More than 50% 
of them have been confirmed to be pre-main- 
sequence objects by measuring their Lithium 
abundances with the help of high-resolution spec- 
troscopy (Covino et al. 1997). We find a com- 
panion star frequency that is lower by a factor of 
0.61 ± 0.27 compared to solar-type main-sequence 
stars (Duquennoy & Mayor 1991). The separa- 
tion distributions are remarkably different: We 
find almost no binaries among the young stars 
with separations larger than 100 AU. 

The binary frequency and separation distribu- 
tion of the confirmed pre-main-sequence stars are 
similar to that of the entire sample, although with 
a smaller statistical significance due to the smaller 
number of systems. 

In contrast to our results for X-ray selected 
stars, Ghez et al. (1997) found a much higher bi- 
nary frequency among 17 stars known before RO- 
SAT. Our interpretation of this discrepancy is that 
the two samples represent two different popula- 
tions. This agrees with the distance determina- 
tions, which indicate that at least some of the X- 
ray-selected stars are only half as distant as the 
stars known before ROSAT. 

It remains unclear how the stars formed. Ac- 
cording to the theory that the multiplicity depends 



on the stellar density, the small number of bina- 
ries and the shape of their separation distribution 
indicate that the stars were once in a much denser 
environment. However, the proper motions show 
no sign of expansion from a single point. This sup- 
ports the hypothesis that the stars formed in small 
cloudlets at their present location (Feigelson 1996; 
Mizuno et al. 1998), although we currently cannot 
explain the bias against wide binaries. 
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